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We propose a new phase of matter, an electronic phase separation transition that starts near 
the upper pseudogap and segregates the holes into high and low density domains. The resulting 
grain boundary potential favors the development of intragrain superconducting amplitudes. The 
zero resistivity transition arises only when the intergrain Josephson coupling Ej is of the order of 
the thermal energy and phase locking among the superconducting grains takes place. We show that 
this approach explains the pseudogap and superconducting phases in a natural way and reproduces 
some recent scanning tunneling microscopy data. 
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The nature of the pseudogap phase has been widely 
recognized to be a key for understanding the physics 
of cuprate superconductors and its complex phase 
diagrampj, (2]. At present there is no consensus on its 
origin and also no agreement on the detailed generic dop- 
ing dependence p of the pseudogap temperature T* (p) [2] . 
This difficulty to find an explanation for the data col- 
lected by many different experiments is certainly due to 
the intricate charge dynamics of cuprate superconduc- 
tors. 

To deal with this complicate charge dynamics we have 
proposed a static phase separation [J13,[E@| based on the 
experimental evidences of ion diffusion in La2CuOi-^-s 
and in Bi2212 above room temperature. The experi- 
mental signals that can be linked with a phase sepa- 
ration are observed at the upper pseudogap T^(jp) (in 
the notation of Ref.fT]) and consequently the ionic seg- 
regation transition must occur at a higher temperature 
Tps(j)). Since T°(j)) ( and Tpsij))) falls to zero in the 
overdoped regime and ionic mobility requires high tem- 
peratures, we assumed previously a charge disorder for 
underdoped compounds and an uniform charge distri- 
bution for p > 0.20[^ 0]- However, new scanning 
tunneling microscopy (STM) data have shown an inho- 
mogeneous local gap structure that remains in the far 
overdoped regimefaH, [13] which cannot be explained by 
an ionic phase separation, due to the low values of T'°(p) 
for large p. 

These STM results on different doping regimes have 
clearly observed local gaps with different amplitudes at 
temperatures below and above Tc{p)^J\^ that ruled out 
ionic phase separation as the sole origin of the cuprate 
inhomogeneities. In order to have an unified description 
of the STM data in the overdoped and underdoped re- 
gions of the phase diagram, we define a distinct phase 
of matter, an electronic phase separation (EPS). In this 
transition, the electrons (or holes) generate bubbles as 
the temperature decreases below the onset temperature 
Tpsip) and freezes at lower temperatures. 

The origin of this novel EPS transition is the proximity 



to the insulator AF phase, common to all cuprates, and 
can be described in terms of competing minimum free en- 
ergy or maximum entropy. As the temperature decreases, 
the entropy of the homogeneous density p becomes lower 
than the anisotropic system made of a granular bimodal 
distribution fsj of AF domains with p{i) « and high 
hole density domains with p{i) w 2p. This condition can 
be written as: 

SM{p)<+Sj^ + SauAp)- (1) 

Where Sm {p) = ^pT is the well known specific entropy of 
a homogeneous fermion gas with density p, Sj^ is the On- 
sager specific entropy for a 2D Ising model with spin cou- 
pling value that yields a Neel temperature at T = 350K, 
taken as a model to the AF phase[ll|. SMix{p) is the 
entropy of mixing p^. In Fig.([T]) we show the condition 
for the EPS transition onset for some selected values of p, 
when the straight lines cross the Sj^ . The used value of 
7 is consistent with the entropy measurements . The 
calculated Tps{p) are in general agreement with the up- 
per pseudogap values^, Q T^{p) and, more importantly, 
it provides a physical interpretation for the origin of the 
electronic inhomogeneities in the cuprates. 

Now that we have discussed why cuprates may go 
through a transition to form granular charge domains, we 
need to describe quantitatively such transition. For this 
purpose we use the theory of Cahn-Hilliard (CII)[12] that 
is appropriate to describe a phase separation transition. 
The difference between the local and the average charge 
density u{i,T) = {p{i,T) — p) is the order parameter. 
Clearly u(i, T) = corresponds to a homogeneous sys- 
tem above Tps{p). Then the typical Ginzburg-Landau 
free energy functional in terms of such order parameter 
near the transition is given by 

/(*, T) = \e^\Vu{i, r)|2 + V{u{z, T)). (2) 

Where the potential V{u,T) = A^{T)u'^ / 4, + 
A^{T) = a{Tps{p) ~ T), a and B are constants that 
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FIG. 1: The lines are S'm(p) — SiHxip) for some p values. 
The intersections with AF entropy yields the onset of EPS, 
that is, Tps{p)- 



lead to lines of constant values of A{T)/B, parallel to 
Tpsip), as shown in the inset of Fig. ([2]). e gives the size 
of the grain boundaries among two low and high den- 
sity phases p±(i)[3, T^. The energy barrier between two 
grains of distinct phases is Eg (T) = (T) / B that is pro- 
portional to {Tps — T)^ near the transition, and becomes 
nearly constant for temperatures close to Tps{p). Thus, 
hereafter we will use Eg{p,T) = V{p,T) as the grain 
boundary potential. V{p,T) — V{p) x V{T) and we as- 
sume, for simplicity, that V{p) have a linear behavior, 
whose equipotentials are parallel to Tps{p). In the inset 
of Fig.© we plot Tps{p), T°{p) both assumed linear and 
the equipotentials are for A{Tps —T)/B. 

For completeness, the CH equation can be written [l5l| 
in the form of a continuity equation of the local free en- 
ergy /, dtu = —V.J, with the current J — MV{5f /5u), 
where M is the mobility or the charge transport coeffi- 
cient. Therefore, 
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A\T)u-B'u'') 



(3) 



We have already made a detailed study of the density 
profile evolution in a 105 x 105 array as function of the 
time steps, up to the stabilization of the local densities, 
for parameters that yield stripe [11 and patchwork 0, [lB| 
patterns. 

The temperature evolution of the second order EPS 
is studied by the ratio A{T)/B^. A{T)/B ^ 0.2 
is close to the value of the measured upper pscudogap 
temperature T^{p) shown in the inset of Fig. ([21). At 
A{T)/B ~ 0.6, the EPS domains are clearly formed as 
displayed in Fig.([2|) and the system is on the limit be- 
tween a disordered metal with grains of two densities, 
and a mixture of metallic and insulator (AF) grains. This 
is possibly the origin of the instability that falls to zero 
near p=0.18 as seen by may experimentsjl, [H, [l3| but 
not detected by the STM data for the Bi2212 seriesfsl]. 
At T w OK the domains are frozen and in general the low 



density insulator regions decrease in number and size as 
p increases, but even overdoped samples have some re- 
maining AF grains according to our simulations, what 
is also experimentally verified by the neutron diffraction 
datajJi]. 
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FIG. 2: (color on line) The charge density map on a 105 x 105 
system after 6400 time steps and A(T,p)/B = 0.6. The inset 
shows estimates of Tps, and some locus of constant values 
of A/B. 

We study the free energy evolution with time and tem- 
perature together with the corresponding density pro- 
file. In Fig. ([21) we show the free energy map associated 
with the density profile of Fig. ([2]), both made by the 
same computer simulation. It shows that the low and 
high density grains, at this temperature (%60 of Tps for 
p = 0.16), are already bound regions of free energy min- 
imum. 




FIG. 3: (color on line) Local free energy (in arbitrary units) 
density profile in the same location and temperature as in 
Fig-©- The dark (red) lines show the potential barrier among 
the grain boundaries. 

As the temperature decrease below Tpg, the poten- 
tial barrier among the grains or the intragrain potential 
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V{p, T) increases and becomes constant at low temper- 
atures. Consequently the holes become confined by this 
effective attraction toward the center of the grains and it 
may be taken as the origin of the superconducting inter- 
action that forms the (intragrain) hole pairs. 

The intragrain superconductivity is naturally study 
with the Bogoliubov-dcGcnnes (BdG) theory in a similar 
fashion as we did before for a phenomenological potential 
and a static phase separation [1, [1, 0]. The calculations 
are performed on a square lattice of 32 x 32 sites, that 
is, on a small part of the charge density profile given in 
Fig.©. 

Assuming the extended Hubbard Hamiltonian to de- 
scribe the hole dynamics, the diagonalization is made by 
the BdG equations [1, 0, with the hopping value 
t = 0.15eV, next neighbor hopping t2 — 0.70i, on-site 
repulsion U = \M and, most importantly, the EPS 
next neighbor attraction V{p, T) derived from the val- 
ues of A{T)/B. Except from the temperature dependent 
V{p,T), all the others parameters are similar to values 
previously usedjl, 0, [1] . 

Following our free energy simulations, from low tem- 
peratures up to Tps{p) when the grains melt down, we 
can obtain the qualitative behavior of V{p^ T). In order 
to yield average coherent gaps values comparable with to 
the STM data on 0.11 <p< 0.19 Bi2212 compounds^, 
we find a set of parameters that can be written as 

Vip, T) = V{p) X V{T) = (-0.9 + 2.8 x p) x 

{l-T/Tps)^^-^'^-^\ (4) 

where the values are in eV ^ ^(p) is linear and vanishes 
at p « 0.32 following Tps{p). V{T) falls to zero near 
Tpsip) and increases towards T ~ OK. 

In general, the CH and BdG combined calculations 
yield very low or almost zero local gaps for the regions 
with low densities, that is, p{i) < 0.09. At the grains with 
larger local densities p{i) > 0.1, the local Fermi level is 
large enough to have d-wave superconducting amplitudes 
Ad(i,T). We define the local superconducting tempera- 
ture Tc{i) as the temperature which Ad{i,T) arises in 
one given site "i". The largest value of Tc{i) in a given 
compound determines the pseudogap temperature T*{p) 
which marks the onset of superconductivity. Since T* (p) 
is close related with the potential V{p, T) it also increases 
in the overdoped region, similar to the Nernst effect [Tgj 
and many other experiments [l|, 

As the temperature decreases below T*(p) and some 
high density grains become superconductors, the zero re- 
sistivity transition takes place when the Josephson cou- 
pling Ej among these grains is sufficiently large to over- 
come thermal fluctuations, that is, Ej{p,T = Tc) ~ 
ksTcip) what leads to phase locking and long range phase 
coherence. Consequently the superconducting transition 
in cuprates occurs in two steps, similar to a supercon- 
ductin g m aterial embedded in a non superconducting 
matrix[20|, first by the appearing of intragrain supercon- 



ductivity and by Josephson coupling with phase lock- 
ing at a lower temperature, what provides a clear in- 
terpretation to the pseudogap phase. Since Tc(p) is not 
directly related with the local or intragrain supercon- 
ductivity, the gaps Ad(«,T) do not change appreciably 
around Tc{p), specially for underdoped compounds that 
have large T* (p) . This fact is verified experimentally by 
temperature dependent tunneling[21] and angle resolved 
photon emissioni22|. 

Using now the theory of granular superconductors [23j . 
Ej{p,T) oc Cn{p) X A{p,T) where Cn{p) is the normal 
conductivity among the grains. As shown in Fig. ([5]) the 
grain boundaries are made of walls with the mean den- 
sity p surrounding the grains. On the other hand, the 
conductivity increases a few orders of magnitude with p 
in the Cu — O plane [241] and Cn{p) is small in the un- 
derdoped region. That is just the opposite average be- 
havior of Ad{i,T) that, following V{p,T), decrease as p 
increases. This gives some insights on the superconduct- 
ing "dome shape" of the resistivity transition with the 
maximum Tc around p = 0.16 in the middle of the EPS 
region {Tps{p « 0.32) = 0). Also Ej oc Jcrf, where Jc 
is the critical current density and is the average size 
of the grains. Taking typical optimum doping values [25|. 
that is, Jc ~ lO^A/m^ and « 50A as one can see di- 
rectly from our Fig.©, we get Ej k. 8meV or Tc « 90K, 
which is a good estimate for the Bi2212 optimum Tc- 

Now we turn to the new STM data that motivated the 
introduction of the EPS concept. We firstly notice that 
the presence of the p « AF insulator {p{i) < 0.03) and 
even low density domains {p{i) < 0.09) which are closer 
to the half filled band and requires a high energy cost to 
accept extra electrons explains why injection of electrons 
produces less STM current than extraction and also why 
this asymmetry increases drastically as p decreases j^. 
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FIG. 4: The BdG calculation for Ad{i,T) at 3 locations on 
the 32 x 32 mesh with average hole doping p ~ 0.24. For 
p{i) < 0.9 the STM signal is from activation over the grain 
boundary potential and Ad{j, T). 

Fig.© shows some of the local BdG calculations on 
selected points over a Cu — O plane as in Fig. (21) to com- 
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pare with the high temperature STM data of overdoped 
{p « 0.22 - 0.24) Bi2212 compounds [1, The smaller 
coherent gaps Ad{i,T) are from different locations on 
metallic grains. The larger gaps originated in the in- 
sulator grains and are due to activation over the grain 
boundary barrier Bij{T) w V{p,T) and the calculated 
Ad{j, T) from a metallic neighbor grain j . As one can see 
in the experimental maps, at temperatures above Tc{p), 
they are always surrounded by a small superconducting 
region d, [13]. 

An interesting consequence of this scenario is 
that the lower density (insulator) grains have larger 
gaps but lower local conductivity that was verified 
experimentally Another consequence is that, despite 
the uncertainty on Tc{i) for very small gaps, the results 
follow close the measured relation 2Ad/ KBTc{i) ~ 80. 



o 













T-45 K 






T=60 K 






T=70 K 






T=90 K 











-0.1 -0.05 0.05 0.1 0.15 

Voltage (V) 



In Fig. ([5]) we use our local gap calculation to plot the 
ratio between the tunneling conductance in the normal 
and in the^uperconducting state as measured by Pasupa- 
thy ct al[lO]. Our calculation is made with their Eq.(2). 
As mentioned above, for the case of insulator grains, the 
total gap is the grain boundary barrier Bij (T) plus the 
Ad{i,T) from a neighbor metallic grain. In Fig.® we 
used \Ad\ = 25meV and Bij{T)=10mev similar to the 
larger gap of Fig.(|4]). Also, in the electronic granular sce- 
nario, due to the charge density oscillations and charge 
tunneling through the grains, the local inverse lifetime of 
the quasiparticle excitations [27| F is an oscillating func- 
tion of the applied bias V and the temperature. As 
shown in Fig.fl]) this phenomenological form of F yields 
the measured density of states [l^l with a structure near 
the applied voltage V = O.leV that is captured by our 
calculations. 



In conclusion we have proposed a new electronic phase 
to cuprate superconductors essentially made of disor- 
dered low and high charge density grains. The grains 
are static at low temperatures but melts slowly and dis- 
appear near Tpsip). Such anomalous phase arises due 
to the proximity of the undoped insulator with AF or- 
der, and may be common to other materials with some 
doping dependent phases, like manganites which possess 
also a pseudogap phase [28|. 



FIG. 5: (color online) The temperature dependent STM signal 
versus applied voltage at a low density grain. The gap 
is due to activation energy over the potential barrier of a 
superconducting metallic grain with a local gap Ad{i,T). 
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